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We study classical solutions in the 5*17(2) Einstein- Yang-Mills-Higgs theory. The spherically 
symmetric ansatze for all fields are given and the equations of motion are derived as a system of 
ordinary differential equations. The asymptotics and the boundary conditions at space origin for 
regular solutions and at event horizon for black hole solutions are studied. Using the shooting 
method, we found numerical solutions to the theory. For regular solutions, we find two new sets 
of asymptotically flat solutions. Each of these sets contains continua of solutions in the parameter 
space spanned by the shooting parameters. The solutions bifurcate along these parameter curves 
and the bifurcation are argued to be due to the internal structure of the model. Both sets of the 
solutions are asymptotically flat but one is exponentially so and the other is so with oscillations. 
For black holes, a new set of boundary conditions is studied and it is found that there also exists 
a continuum of black hole solutions in parameter space and similar bifurcation behavior is also 
present to these solutions. The SU(2) charges of these solutions are found zero and these solutions 
are proven to be unstable. 
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I. INTRODUCTION 



The first regular solution of the Einstein- Yang-Mills (EYM) theory found by Bartnik and McKinnon 
[jj two decades ago stimulated intensive research into Einstein-non-Abelian gauge theories because of 
their rich geometrical and physical structure. Soon after, the black hole solution and the violation of the 
no-hair conjecture for these black holes were discovered 0, 3|. Because of this, this model was extended 

Einstein-non- 



5 L Einstein- Yang- Mills-dilaton 
9j theories (see 10| for review). Black holes in 



to a few other models, such as a Einstein-Skyrme ^ 
Abelian-Proca 0, Q and Einstein- Yang-Mills Higgs 8, 
these theories also violate the no-hair conjecture, and some of them are even claimed to be stable. 

Among these, it is found in Ref. [§} that in the spontaneously broken phase of the EYMH model 
the black hole solution can possess a non-trivial field structure outside horizon. However it is not known 
whether there exist other solutions (in particular, black hole solutions that violate the no-hair conjecture) 
to this EYMH theory and whether these solutions are stable or not. In this paper, we extend the work of 
Ref. Q by thoroughly studying other non-magnetically charged and asymptotically flat solutions of the 
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EYMH theory with a Higgs doublet, a Higgs potential and a cosmological constant. Besides the solutions 
in Ref . Q| , we also found two other set of regular solutions and one more set of black hole solutions in this 
theory. One set of the regular solutions and the black hole solutions are for the minimal EYMH model 
(with scalar potential equal to zero) and the other set of regular solutions has an oscillatory decaying 
feature. The black hole solutions are characterized by the nodes of fields and therefore also provide a 
counter example for the no-hair conjecture. All these solutions are proven to be unstable. Therefore if 
one demands stability in the no-hair conjecture, the black hole solutions found here will not conflict with 
it. 

The paper is organized as follows. In section [Til we first describe the EYMH model. Then we give the 
spherically symmetric ansatze for the metric, Yang-Mills field and Higgs field, and derive the equations of 
motion. In section iHll we study the boundary conditions and asymptotic behavior that are compatible 
with the field equations and certain physical criteria. In sections IIVI and [Vj numerical solutions for 
both regular and black hole boundary conditions are found using the shooting method, and their general 
features are discussed. We emphasis that all solutions found here are new solutions that were not found 
in previous studies. In the last section IVT1 we discuss the stability of the solutions and discuss possible 
extensions. 



II. MODEL, EQUATIONS OF MOTION AND ANSATZE 

The action of the EYMH model could be written as 

S = J d 4 x^(£ g + £ m ) (1) 
where g = det g^ v and the Lagrangian densities of the gravity part C g and the matter part C m are 

£ * = I^ + A ' (2) 

C m = -\F$FQg*<f"> -[D^[D^]g^ -V(&$), (3) 

T/($t$) = m 2 $t$ + A ($t$)2 ( 4 ) 

The vector field is of SU (2) and the doublet Higgs field <J> has two complex components: 



where r a /2 (a = 1,2,3) are the SU(2) generators and (j>i and </>2 are complex fields. The field strength 
Fftu = d^Ay — d v A^ — ig[A^ 1 A v \ is defined as in flat spacetime and = — igA^ is the covariant 
derivative, where g is the gauge coupling constant. R^ v and R are the Ricci tensor and scalar respectively, 



and m and A are the mass and coupling constant of the Higgs field. A appears as a parameter to assure 
the positive definiteness of the energy. 

Variation of the action ((T|) with respect to g^" , Av and $^ gives rise to the following equations of 
motion 



^Gn U$F(fg»> - \ gia ,Fj$F<#g*°gn + 2(D^(D U $) 
d^(V^gF^) + J—g~ge abc A^F^ - 2Tr[(d> i T a )T b }F$ 



Df, [y/^jg^Dv®] - 7^[m 2 $ + 2A($ t $)$] = 0. 



fJLV 

= o, 



(6a) 

(6b) 
(6c) 



The first equation could be written in the trace-reversed form 



— RitGm < F^ F^o Xp — -a F^ F^ a Xa a m 
v — 07h_tat <r x r v g g pu r x r g g 



(7) 



To find useful solutions to this equation system, we concentrate on the static and spherically symmetric 
ansatze of the metric, gauge field and Higgs field. The most general form for static and spherically 
symmetric metric could be written as follows 



ds 2 = -T{r)-' z dt z + R{r) 2 dr 2 + r 2 (d9 2 + sin 2 6 dtp 2 ) 



(8) 



where if we let 



R(r) = (1 -2M{r)/r)- 1 / 2 1 



(9) 



M(r) could be interpreted as the Misner- Sharp mass within radius r [12|. 

For the study of black hole solutions, it is more convenient to use an alternative metric of the form 



ds 2 



(l - ^M) e - 2 *(^ 2 + (l - _1 d, 2 + r 2 (d6 2 + sin 2 9dip 2 ), 



(10) 



which is obtained from Eq. ((5|) by defining 8(r) = — ln(i?(r) /T(r)). A regular event horizon of this metric 
requires that 



M{r h ) = -f, S(r h ) < oo, 



(11) 
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where is the horizon radius. For later purpose, we can always rescale the time coordinate in (JSJ and 
(fTOj) by dt = T- 1 (0)dt and dt = e~ 5 ^dt respectively, after which we have f(0) = 1 and S(r h ) = 0. 
These rescalings will simplify the numerical calculations that will be done in the next section. Hereafter 
we will drop the tilde symbols. 

The most general spherically symmetric form of gauge field is given by Q, Q, [l3 1 



A = -z{a,T r dt + br r dr + [dr e - (1 + c)T v ]d6 + [(1 + c)t 8 + dr v ] sm%}, 



(12) 



where a, b, c, d are real functions that only depend on r and t, and (r r , tq, tJ) are the Lie algebra su{2) 
bases which satisfy tr(T T&) = l/26 a b, [ T a,i~b] = ie abc T c (a,b,c = r,6,<p). In particular, here we choose 
them to be 
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~ 2 
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{r 3 J 







(13) 



where J is just the unitary Jacobian matrix for transformation from cartesian to spherical coordinate and 
tj (i = 1, 2, 3) are the usual Pauli matrices. The gauge field 1)120 has a residual £/(l) gauge transformation 
of the full SU(2) gauge group 



U = exp(i(3(r,t)T r ) 

A -> UAU- 1 + -UdU~\ 
9 



(14) 
(15) 



where j3(r, t) is an arbitrary function. This gauge transformation could be used to put b = in Eq. JT! 
identically. For the remaining three degrees of freedom, we eliminate two by concentrating only on the 
purely magnetic and static YM field, i.e., setting a = d = and c = c(r). With this setting, the ansatz 
(fT2l is reduced to 



A = - [-(1 + c)r v d9 + (1 + c)t 9 sin9dip] . 



The most general form of the Higgs field could be written as 



(16) 



s 1 / ^2 (a) + itpx(x) 
V2 I ^(a;) - # 3 (a;) 



(17) 



where we can treat the three degrees of freedom ip a (a = 1,2,3) as those of a vector field tp(x). To 
yield a spherically symmetric and static energy density, we can use a more useful ansatz by taking 
4>{x) = 4>{r), tj}{x) = ip(r)h r To simplify the numerical analysis that will be conducted in the 
following sections, we further set ip = henceforth. 
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To obtain the detailed equations in terms of the fields in these ansatze, one can substitute the ansatze 
{5}, p^|) and <j>(x) into the system ([5]). It is found that the system of equations of motion consists of 
non-linear ordinary differential equations of w (w(r) = c(r)), <p and M 

r (l - —) w " = ( (m 2 2 + - 2A)r 2 + ll + wfcj? + 1(1 - w 2 ) 2 - —) w' 

\ r J \ 2 2 r A r J 

1 . , , w(w 2 — 1) . 

4 r 

1 _ V = f (m 2 2 + 1 A0 4 _ 2A)r 2 + 1 (1 + w)2(f) 2 + 1 (1 _ ^2)2 + 2M _ 2 



r / \ 2 2 r 2 r 



_ A r0 3 + + mVl 0, (18b) 

/ 2M\ / , 2 1 2 \ 1(1-^ 1 
M = J (™ +2^) J + 2 + 4<M 1 +™) 

imV + W- A V> ( 18c ) 



2 ^ 4 
and for the regular solutions, of T 



2M\T> ( 2M\ ( a 1 2 \ l(L-u?f 1 2 2 



r / T V r y \ 2 v r ' / 2 r 2 

1 2,2 , 1, , 4 A \ 2 M 



_ m - r + _ A ^_Ajr 2 -— , (19) 

where ' refers to derivatives with respect to r. Note in these equations and hereafter, we set g = 1 and 
Gat = l/(47r), which are equivalent to the rescaling r — > r-y/G/y/g and 5 — >■ gv / G/v"S', and therefore the 
system on longer depends on them. 

For black hole metric (TO]), we replace T(r) in Eq. (JTScJ by T(r) = e*( r >(l - 2M{r)/r) 1 / 2 and Eq. ([II]) 
for T(r) is replaced by the field equation for 8{r): 

2m' 2 

8' = -rcb' 2 - —. (20) 



Note that the Eqs. (fT8"] ) -([20 l) are identical to equations (4.19)-(4.24) in Ref. [8| and equations (17)-(22) in 



Ref. 



ll|. However here we will solve these equations by setting parameters differently. If <p(r) is put to 



zero, then equations (3)-(5) in Bartnik and Mckinnon's work [jj could be recovered. Because the system 
has mirror symmetry 4>{r) — > — </>(r), we only need to concentrate on one case (we pick <fi(r — > 0) < 0) 
for the specification of the boundary conditions, which are necessary in order to solve the equations 
numerically. 
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III. BOUNDARY CONDITIONS 



The system of equations (|T5|) and (fPT))) and the system (JTBJ) and (f2T))) are only solvable numerically 
provided proper boundary conditions are given. The former system has singular points at r = and 
r = oo and the latter has singular points at r = r/, and r = oo. In the following two subsections we 
discuss these boundary conditions. We first examine the behavior of regular solutions at the singular 
points. 



A. Boundary conditions for regular solutions 



As r — > oo, we require that the gauge field w(r) and scalar field <f>(r) approach constant asymptotic 
values while the metric function T(r) goes to a constant but nonzero value. 

For asymptotically flat spacetime, we can do a careful asymptotic analysis by substituting the following 
expansions 



w(r) 
0(r) 
M(r) 
T(r) 



— 1 + 6w(r), 
4>a + 5</>(r), 
M + SM(r), 
T + 8T(r), 



(21a) 
(21b) 
(21c) 
(21d) 



into the Eqs. (jT5J) and (|19p and treating Sw, 6<fr, 6M, 8T as perturbations. Here we know that the 
asymptotic value for w(r) is -1 and <f>o, Mq and To refer to some asymptotic constants. 
It is found that there are three possible asymptotics. The first takes the form 



w{r) 1 + ce*" 72 , (22a) 

1 cV 01 " 



W~^o + ^— 5-, (22b) 



M(r) ~ M + ic 2 o e 0or , (22c) 

T(r) ~ To + ^p, (22d) 

where c is some positive constant, with the following conditions on the parameters in the model 

A = 0, ^ 0, m 2 = A = 0. (23) 



The above parameter setting means that the scalar potential is zero and there is only a kinetic energy of 
the scalar field that is coupled to gravity. Therefore this case corresponds to the minimal EYMH theory. 



7 



The second is given by 1 



c 



w(r) 1 + -, (24a) 



r 

4>(r) ~ — — (ci sin (\/ — m 2 r^ + c 2 cos (\/ — m 2 r^ ^ , (24b) 

M (r) - M + y dr-^ cos ^V^m^r) , (24c) 

r (r) ^T + ^, (24d) 
r 

where c and ci, C2 are some constants and p(r — > oo) > 1 is a slowly increasing function (slower than 
r 1 ), with the condition A = 0. 

The last is the asymptotics that satisfy 

<ft = =! f, A = -IX' (25) 

The equation system with the last asymptotics has been solved in Ref. Q|. Therefore we will only 
concentrate on the first two cases. We emphasis that the above three cases are all the allowed asymptotics 
for the EYMH theory under the demand that the solutions are asymptotically flat. 

For r = 0, regularity of the metric and fields requires that there exists a series solution to each function 
and in particular T(0) ^ 0. We also impose the finiteness condition at r = for the energy density Too, 
which could be calculated from the right side of Eq. ([7]). Using the field equation (TT5|) . one can show 
that Too satisfies 

M'(r) 

T oo(r -4 0) ~ 5^ (26) 

and thus the finiteness of Too(0) requires M(r) ~ 0(r 3 ) as r — > 0. These conditions directly lead to the 
following two possible sets of boundary conditions at small r, classified according to the number of nodes, 
denoted by k (k G Z + ), of w(r). For odd-fc solution, we have 

w(r) = 1 + ar 2 + w 4 r 4 + 0(r 6 ) (27a) 

<j>(r) - hr + for 3 + 0{r 5 ) (27b) 

M(r) = M 3 r 3 + 0{r 5 ) (27c) 

T(r) = f + T 2 r 2 + 0{r 4 ), (27d) 



1 The function p(r) in I l24bl and H24cl l is given by a complicated differential equation that can not be solved analytically. A 
numerical inspection indicates that it is a function that increases slower than r 1 . The asymptotics for M(r) is obtained 
by setting c\ = for (j>(r). The general form for M(r —> oo) without setting ci (or C2) to zero is given by an equation 
that dose not allow us to get a compact solution. However, the oscillatory feature for M(r) in this case should still be 
present. 
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where 



f/>4 



— (16a 3 + 6a 2 + 8a(6 2 - A) + b\) 



•>z = 



10 



(36? 



2 , 2 

1 m 



1 



T 2 = -2a 2 - -A 



2a) 



15 



-A 



(27e) 
(27f) 

(27g) 
(27h) 



For even fc solutions, the boundary condition is 



w(r) = -1 + ar 2 + w 4 r 4 + C(r 6 ) 

0(r) = & o + 2 r 2 + O(r 4 ) 

M(r) = M 3 r 3 + 0(r 5 ) 

T(r) = 1 + T 2 r 2 + 0(r 4 ), 



(28a) 
(28b) 
(28c) 
(28d) 



where 



M'4 
<l>2 

M 3 
T 2 



h 3 ~ To a2 + h (1 + 4A6 ° + 8m2) 6 ° a " l Aa 

1 



(A6 2 + m 2 ) b 



2a 2 



-2a' 



12 



2m 2 



12 



2m' 



Mg) - 3 A 
+ A6 2 ) b 2 - X -K. 



(28e) 
(28f) 
(28g) 
(28h) 



B. Boundary conditions for black hole solutions 

For black hole solutions, we have to set 

M(r h ) = r ±, 5{r h ) = (29) 

and require the gauge and Higgs field to be finite and smooth at r = r% in order to have a regular event 
horizon. Using Taylor expansions for fields w(r), <f>(r) and M(r) near the horizon, we can get a valid set 
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of boundary data at r = : 



r h - (1 - w 2 h ) 2 /r h - 2 (4> h /2) 2 (1 + w h yr h - (m^ , + A0 4 /2 - 2A)rg 

x = (l + W/ ,) 2 ^/2 + (m 2 + A^)^rg 

r ft - (1 - t^) 2 /r fc - 2 (<K/2) 2 (1 + w h )*r h - (m 2 ^ 2 + A0 4 /2 _ 2A)r 3 ' 



M'(r h ) = \ (m 2 4>l + \Mh - 2A) rl + 1(1 + w h ) 2 ^ h + 1(1 - w h f/r 



2 / 2 
ft) 



(30b) 
(30c) 



where w{rh) = and <f>(rh) = 4>h- At r = oo, the same asymptotic analysis as in the previous subsection 
can be done. It is found that for the case (f2"2"T) . 8{r) takes the form 

*(r)~<y -i^— , (31) 
2 r 

and for the case (f2"4")l 

S(r)~5 -Jdr r(^(r)') 2 , (32) 

where <$</>(r) is given in (|24b|) . 



IV. REGULAR SOLUTIONS 

With conditions (|2T|) or (|28p at r = 0, we use the shooting method to solve the Eqs. (fT5|) and (fT9")l 
numerically. Using a standard ordinary differential equation solver, we evaluate the initial data for the 
functions at r = 10~ 2 and use tolerance 10~ 12 to shoot the parameters (a, bi) for conditions (f2~T)) and 
(a, bo) for conditions (f2"5)) and integrate towards r = oo to match asymptotics (|22l)(l2"4")l . A drawback of 
the two parameter shooting procedure is its slow convergence. Therefore, we will limit our study to only 
k = 1 and/or k = 2 solutions. Scalar mass m, scalar coupling A and cosmological constant A are the 
three parameters on which the shooting process depends. We will clarify how they affect the existence 
and features of the solution. In the following, subsection IIV Al contains solutions that match asymptotics 
(|2"2"|) . while in subsection IIV B[ solutions match condition (pH)) are shown. 

A. Asymptotically flat solutions with asymptotics (122|) 

For system with both boundary conditions (|27p and ([2"5]) at r = and (|2"2"|) at r = oo, we found that 
when &i,o(&i or bo)=0, there exist only solutions with a taking discrete but infinitely many values whose 
magnitude falls in the range of 0.453 and 0.707. We denote these values by a*, where k are positive 
integers that equal the numbers of nodes of w(r). Indeed these solutions are just the solutions of EYM 
theory discovered in Ref. [l| and a list of values of and the position of the nodes of w(r) can be 
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FIG. 1: The k = 1 and k = 2 parameter curves. Each point on these curves corresponds to a k node solution. 
The circle(s) at 61 = is the a\ and at bo = are ai and 0,2 of EYM theory. 

found in Ref. [jj. In these solutions, 4>{r) is identically zero and this identically vanishing <fi(r) is a good 
solution because it dose not contribute any kinetic or potential energy to the system and therefore allows 
the possibility of asymptotic flatness. 

The solutions of our EYMH theory ([1]) and that are unique to boundary conditions (|2"T)) or (|2"5)l are 
those with bi^ 7^ 0. For each positive integer k, we found that there exists in the shooting parameter 
space spanned by a and b\ ,0 a continuous curve, each point of which can give a k node solution. The 
two end points of this curve are the (a*,, bi t o = 0) and (a&_i, 61.0 = 0), while the points in between have 
afc < a < afc_i and &i,o 7^ (for k = 1, 2 curves, see Fig. QJ. 

The main feature of these solutions is that there is a bifurcation of the solutions, with bifurcation points 
being the peak points (points P in Fig. [1} of the parameter curves. On each parameter curve, from the 
peaks along the negative direction of \a\ to point (ak-i, &1.0 = 0), one find that the configuration of the 
solutions resembles that of the k — 1 node EYM solution; from the peaks to (a/,, 61,0 = 0), the solutions 
resemble the k node EYM solution (see Fig. 0); while the solutions at the peaks P have the feature 
that their asymptotic mass (or energy) is relatively larger than other solutions of the same parameter 
curve. This variation of the solutions along the parameter curve manifestly shows that the solutions to 
the EYMH theory naturally cover that of the EYM solution, as one should expect. 

This bifurcation is very different both superficially and in its origin from the bifurcation behavior 
discovered in Ref. There the bifurcation depends on the variation of one original parameter in the 
theory - the mass of scalar, while here the bifurcation occurs purely in the shooting parameter space. For 
the origin of the phenomenon, in Ref. Q| the EYHM theory is not minimally coupled: they have scalar 
mass and (j) 4 potential terms, and the bifurcation was explained using two length scales L\ = g~ 1 G 1 / 2 



and L2 = g^ 1 \J X/m 2 . However here in this subsection, the Higgs field is minimally coupled to the gauge 
field and gravity and therefore essentially there is only one length scale L\. Therefore we can not use the 
same length scale argument and we tend to believe that this bifurcation is simply due to the variation of 
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Higgs field <f>(r) and gauge field w(r) along each parameter curve, which is purely an intrinsic structural 
feature of the EYMH theory under consideration. 

This is supported by the observation from equation (|18cp (with fj, = m = A = A = 0) that the mass or 
energy of the system is crucially dependant on the derivatives w(r)' , <fi(r)' and the deviation of w(r) from 
±1 and <f>(r) from 0. The derivatives represent the kinetic energy from the gauge and scalar field; and the 
other terms come from the covariant derivatives of the scalar field, where the non-Abelian field couples 
to itself and the scalar field. (We emphasis that this feature is unique to non-Abelian gauge theories.) To 
see how the parameters a and 61,0 determine the solutions in Fig. [3J we take k = 1 for example. Along 
the directions of the parameter curve that |6| increase, from (|27h|) one see 0(0)' increases while w(0)' dose 
not change. This effectively increase the rate that M(r) grows at small r and eventually leads to larger 
asymptotic mass at the peak P of the parameter curve. The process that the fields and mass function 
evolve to larger values of r has to be determined from the field equations (|18 [l -(|19 [l . whose complicated 
structure blocks us from gaining further insights about the physics. 

Note that throughout this subsection, the solutions found are the solutions to the simplest EYMH 



12 



models in the sense that both scalar-gravity and scalar non-Abelian gauge field are minimally coupled. 



B. Asymptotically flat solutions with asymptotics (1241) 



From the asymptotics (|24bjl and (|24cj) . it is seen that that the asymptotic solutions to <j>(r) and M(r) 
will be oscillatory with a decreasing magnitude, provided that m 2 < 0. The situation for m 2 > will 
produce hyperbolic asymptotics (from (|24b|) and (|24c[) h which corresponds to a spacetime that is far 
from being flat and therefore is not of interest here. Therefore in this subsection, we will always study 
the case with m 2 < 0. In solving the equation system with asymptotics (|24[) . we will simply set A = | in 
order to reduce the amount of calculation. 
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FIG. 3: The solutions for k = 1 to the system (|18[) and (|19[) with asymptotics (J24J) . For all these solutions, 
A = 1/8, fj, = and A = 0. For the lower left solution, m 2 = —0.04. For all other solutions, m 2 = —0.01. 



We study the solutions satisfying (j27|) with k = 1 first. When in 2 is less than zero, it is found that 
two continuous curves in the parameter space spanned by a and b\ start to emerge form the point 
(a = ai,bi = 0) and (a = 0, b\ = 0). These curves are extended as m 2 decreases and finally join each 
other at m 2 w —0.046. Each point on these curves can give a valid solution that satisfies the boundary 
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FIG. 4: The k = 1 parameter curves of the system ()18p and (|19p with asymptotics (|24p for different m 2 's. The 
two lower end points of the left figure are (a = -0.4228, bi = -0.0975) and (o = -0.02438, bi = -0.03467) and 
center figure (a = -0.2944, 6i = -0.1914) and (a = -0.1288, 6i = -0.1474). 




FIG. 5: The parameter curve and solutions for k = 2 to the system with asymptotics (|24p . For all these solutions, 
m 2 = -0.002, A = 1/8, /i = and A = 0. 



conditions. These solutions in general are oscillatory and only approach their expectation values at 
infinity However, as we decrease b\ (increase |6i|; see Fig. 0|) along these two solution curves to their 
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end points, the starting radius of the oscillation of 4>{r) becomes delayed and the first minimum of <f>(r) 
is enlarged, and eventually at the end point of each parameter curve, 4>{r) becomes flat right after it 
reached its nonzero minimum value. This violates the asymptotics (f2"4"]l and will force M(r) to diverge to 
negative infinity and therefore the end points do not correspond to physical solutions (see Fig. 0] lower 
right). For M(r), with b\ decreases (|6i| increases), its oscillation amplitude grows larger until b\ reaches 
its end point value, where the oscillation amplitude of M(r) becomes infinity. Fig. [3] shows the solutions 
for m 2 = —0.01 and various values (a, &i)'s. Fig. 0] shows the extension and joining of the parameter 
curves with respect to the decrease of m 2 . 

The phenomena where the parameter curves are extended with larger |m 2 | and where the solutions are 
oscillatory could be understood in a heuristic manner. The key is still the field equation (|18cj) of the mass 
(or energy) function. The existence of asymptotically flat solutions for M(r) depends on the balance of 
the kinetic term ^(rcj)') 2 and mass potential term |m 2 ^ 2 . As m 2 (< 0) decreases, the mass term provides 
a larger negative value so that the allowed <j>' could also be extended to a larger value, which means b\ 
is extended noticing Equ. (|27b|) . While for the oscillation of M(r) and the increase of the oscillation 
amplitude, it purely comes form the mass term ^m 2 (f> 2 because only this term is negative in M{r)' and 
larger |m 2 | provides larger slope for oscillation of M(r). These are confirmed by the simultaneity of 
oscillations of M(r) and 4>(r) in the solution shown in Fig. [3J The explanation of the divergence of M(r) 
for solutions corresponding to the lower ends of the parameter curves (left and center of Fig. U) relies 
not only on Eq. (|18c|) but (|18b[) . From the latter we see that if 4>(r) becomes flat at some large r, one 
has to have Atfi 2 = — m 2 , that is, <f> = \J—m 2 /\ and this means that there is a negative energy in (|18c[) 
M(r >> 1)' = — m 4 /(4A)r 2 , which directly lead to the divergence of M observed in Fig. |3j (lower right). 
The bifurcation of the solutions for each m 2 < —0.046 is similar to that of the solutions discussed in the 
previous subsection. Again, because the bifurcation is for different choice of the inner parameters a and 
6o.i but not with respect to different choice of m 2 , we can not explain the bifurcation by appealing to 
the existence of the two scales, as was done in Q|. 

For solutions with even nodes (see Fig. [5]), we similarly found that when > m 2 > —0.002, there exist 
two curves in parameter space spanned by a and bo and eventually these two curves join at m 2 w —0.002. 
These solutions also have oscillatory <j){r) and M(r) and their dependence on various parameters are quite 
similar as that of k = 1 node solutions. 

V. BLACK HOLE SOLUTION 

Here we solve numerically for the boundary conditions (|2"2"j) and ([30)) using Wh and </>/j as shooting 
parameters and letting r/j = 1. We evaluate the boundary data at r = + 10~ 2 first and then integrate 
using 10 -12 as tolerance towards r — oo to match asymptotics (|2"2"|) . The solutions to these conditions 
only have kinetic energy contribution to the total mass from the Higgs field. Again, the solutions exist and 
bifurcate along continuous curves in the parameter space spanned by Wh and <ph (see Fig. [5] for parameter 
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FIG. 6: The parameter curves for k — 1 and k = 2 solutions. The blue circles are the Wh(l) and Wh(2) of EYM 
black hole solutions. 




FIG. 7: The k — 1 and k — 2 black hole solutions. The solid curves are solutions corresponding to peaks P in the 
parameter curves. 
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curves and Fig. [7] for solutions). The end points of these curves form a discrete set {(u)h{k), 4>h = 0)} 
where k is the number of nodes of w(r). The solutions corresponding to these end points indeed are 
just the black hole solutions of EYM theory found in Ref. Q and a list of the first five Wh{k) could 
also be found there. For the points (wh,4>h) on the curve that satisfy Wh{k) < Wh < Wh{k — 1), the 
corresponding solutions of w{r) have k nodes. Similar to the regular solutions, the black hole solutions 
for each k in general also bifurcate into two classes: the k node solutions that resemble the k — 1 node 
black hole solutions of EYM theory as \wh\ increases along the parameter curve and the k node solutions 
that resemble the k node black hole solutions of EYM theory as \u>h\ decreases. Again we take the 
bifurcation point to be the peaks of the parameter curves, where the solutions have asymptotic masses 
that are relatively larger than solutions with smaller \<f>h\- This bifurcation behavior is again believed 
to be due to the inner structure of the field equations; but unlike the regular solutions case where we 
have a simple relations (f2~T|) and (|2"5|) between the shooting parameters a, bi o and <fi'(0) and u>'(0), the 
relation (|30p between Wh, 4>h and w'(rh), 4>'{rh) are quite complicated. We therefore would not study 
this bifurcation in more details, but just to remind the readers that the solutions here are the black hole 
solutions to the EYMH theory with scalar field minimally coupled to gauge field and gravity. 



VI. DISCUSSION 



An important global property of these solutions is their charge. The non-Abelian su(2) electrical charge 
Qe and magnetic charge Qm of the gauge fields can be defined as [3| 

7Z I dS ^V-g | *\ 



Qm 1 

Using the given ansatze and the asymptotic values, it is found that for all type of solutions 

Qe = 0, Q M oc (1 - w{^) 2 )T r = 0, (34) 




which means the solutions are chargeless in the gauge corresponding to ansatz ([T| 

Another important issue is the stability of these solutions. For the equation systems with conditions 
(|25p . we know that both their regular and black hole solutions are unstable with respect to linear per- 
turbations [15I. fl^. Because the stability depends crucially on the boundary conditions and asymptotics, 
we need to do a separate examination of the stabilities of the solutions found in this paper. We can carry 



out t he p erturbation of SU (2) gauge fields in fT2")) and Higgs field in (fTT)) along the line of Ref. 15| and 
Ref. |l6| (note we have a sign difference in the definition of w(r) compared to Refs. Q, H|): 



b(r,t) = -w'(r)z(r)e iut , d(r,t) = [w(r) 2 - l]z(r)e luJ \ i/)(r,t) = -0(r)z(r)e lut . (35) 
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For regular solutions, we set z(r) = 1. For black hole solutions, we choose z(r) to be a real function that 
will be determined later. 

Using this perturbation, it is found that the perturbation equation take the same form 



(36) 



where ^ = (b(r, t),d(r, t), -0(r, t)) T and H and A are matrix operators as in Refs. [I5l.ll6|. The expression 
for the eigenvalue square 



is also the same as in Refs. 
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3 



(37) 



mm = 
(mm = - 

+ 



r(w') 



'12 



ro 
oo 



r 
oo 



(w 2 -l) 2 (w+ l) 2 2 r 2 "l 



S 

2N{w'f 



2N{w') 2 + 2 



NS 

{w 2 - l) 2 



4NS 



z{r) 2 dr 
Sdr 



(38) 



(1 - z{r) 2 ) Sdr 



2(w 2 



l) 2 + |(» + lW 



(z') SNdr 



2(w 2 



i) 2 + J» + i)W 



(39) 



where N = 1 — 2M/r and S* = (1 — 2AI/r)~ 1 T~ 1 , ro = for regular solutions and 5* ~ e , ro = rh for 
black hole solutions. The terms in the last line of Q39p are the boundary terms that were only impli citly 
mentioned in Ref. [l^. Because the expressions (|3"5)) and ([59^1 are exactly the same as in Refs. 15 



our arguments will follow exactly these references. We only need to pay attention to the applicability of 
our asymptotics ([2"2"| , (|2"4"|) and |3T|) and (|3^|) , because they are the only relevant difference from those of 
Refs. QQ. Below, we show the details of the arguments. 

For regular solutions (z(r) = 1), it is clear that the last three lines in (|30|) drop and therefore (4'|i7|il') 
is clearly negative definite. From the asymptotics (|2"2"|) and (|2"4"|) . with little effort one can see that as r 
becomes large all terms in the integrand of (\E r |A|\t r ) vanish at least at the speed of 1/r 2 and therefore 
(SI/ 1 ^4 1 \E f ) is finite. This shows that the eigenvalue is negative and therefore the solutions are unstable. 

For black hole solutions, we will show that for some properly chosen z(r), (*E f |i/|\t r ) is negative definite 
and terms in the integrand of 1 ^4 1 NI/ ) decrease fast enough so that this term is finite. Moreover, we also 
need z(r = rh) — for the perturbation (| 35[) to vanish at the horizon. These conditions can be satisfied 
by a z(r) chosen according to Ref. ly, |l_7|. This involves defining the tortoise coordinate r* by 



dr* 
dr 



1 

~NS 



(40) 
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and a set of functions Zk (r* ) by 



( 



(41) 



where 



-D < 




< for r* e [a,a+ 1], 



1 for r* S [0,a], 



u(r*) 



for r* G [a+ l,oo]. 



(42) 



If we let z(r) = Zk (r* (r)) in (f3U)) , we see that the boundary terms (last line) vanish at horizon (corresponds 
to r* = — oo) because Zk(r* = oo) = and vanish at space infinity because of the asymptotics (|2"2"j) . (|2"4")l 
and pi]) and (f3"2"1) . One can also see that the second and third lines of will vanish uniformly as 
k —> oo. This establishes that (\J r |.H'|\f') is negative definite. This choice of z(r) with \z(r)\ < 1 and the 
asymptotics ((22|) . (|24p and (|3Tjl and (|32| also guarantee that the integrand of ([38)) vanishes at least at 
the speed of 1/r 2 as r becomes large and therefore (^1^1^) is finite. Therefore again the eigenvalue is 
negative and the solutions are unstable. 

Even though the EYM theory and its extensions have very interesting theoretical features, they are only 
relevant in few known objects in nature. One of these is the neutron stars, where the matter becomes 
very dense so that the gravity is very strong and the gluon fields become the fundamental degrees of 
freedom. The matter in neutron stars is usually described by effective Quantum Chromodynamical 
theories and chemical potential in this matter is large and important. Our original hope for this paper 
was to study the effect of a nonzero chemical potential to the EYMH theory. However, the reduced 
spherical symmetric form of the gauge fields - ansatz (|16p - that is used by previous studies on EYMH 
theory, dose not allow a consistent and simple introduction of the chemical potential that respects the 
spherical symmetry. Because a non-spherically symmetric setup in this case is very difficult, one can only 
attempt to construct a consistent theory with nonzero chemical potential and spherical symmetry from 
other forms of reduction of the most general gauge field ansatz ([T2|) and its gauge equivalents (e.g., see 
Ref. [3|). The result of this attempt will have to be reported later. 
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